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ABSTRACT 


Beginning  with  Boltzmann's  equation  and  the  equation  of 
momentum  transfer  for  a  partially-ionized  gas,  seven  basic 
equations  are  derived  for  small -amplitude  perturbations.  These 
are  used  to  obtain  a  determinant  whose  eigenvalues  identify  five 
modes  of  propagation  in  the  most  general  case.  A  steady  mag¬ 
netic  field  and  an  arbitrary  direction  of  propagation  with  respect 
to  this  field  is  supposed,  and  the  complications  introduced  by 
this  general  approach  are  examined  by  a  brief  study  of  simpler 
cases. 

It  is  shown  that  only  by  means  of  a  high-speed  electronic 
computer  may  attenuation  factors  and  phase  and  group  velocities 
be  accurately  determined.  This  paper  surveys  work  done  by 
other  investigators,  includes  the  approach  taken  by  a  research 
group  at  the  Electrical  Engineering  Research  Laboratory  which 
has  access  to  the  Control  Data  Corporation  1604  Computer,  and 
is  only  to  be  considered  preliminary  to  later  presentation  of 


computer  results. 


I.  INTRODUCTION 


The  dispersion  relation  is  an  equation  which  can  be  used  to  calculate 
phase  and  group  velocities  and  attenuation  and  phase  factors  in  a  medium  of 
propagation.  It  is  well  known  that  many  parameters  enter  into  this  equation. 
The  presence  of  a  fixed  magnetic  field,  a  situation  existing  in  the  ionosphere, 
for  example,  greatly  affects  the  form  of  the  dispersion  relation.  This  paper 
will  outline  the  methods  used  by  various  investigators  in  an  attempt  to  arrive 
at  a  general  dispersion  relation,  without  making  the  various  approximations 
which  usually  make  the  equations  more  tractable. 

Three  types  of  waves  have  been  identified  in  ionized  gases  - 
electromagnetic,  hydromagnetic  and  electro-acoustic  or  space  charge  waves. 
Dispersion  relations  for  the  fully-ionized  gas  have  been  obtained  under  the 
special  assumptions  which  make  possible  writing  down  neat,  separate 
relations  for  each  of  the  three  types.  However  it  has  long  been  desired 
to  find  a  relation  which  is  most  general  and  reduces  to  the  simpler  relations 
for  each  case. 

We  pose  the  problem  as  follows;  Axes  are  set  up  in  an  ionized  gas 
of  three  components  -  electrons,  ions  and  neutral  particles.  We  suppose  a 
fixed  magnetic  field  pointing  in  the  positive  z  direction,  B  =  k  A  plane 

wave  front  is  supposed  as  shown,  0  being  the  angle  between  the  propagation 
direction  and  the  B  vector.  Consider  all 


-  I  - 


2 


+x 


Figure  1  Co-ordinate  Axes 

derivatives  with  respect  to  the  y-axis  equal  to  zero.  The  quantity  "t)"  in 
the  e  ^  term  is  the  unit  distance  along  the  propagation  direction.  A  field 
quantity  is  then  written 

A  =  A  e"^  e^"*  =  A  sin  0  +  z  cos  9)  ^jwt 

o  o  ' 

j  wt 

and  hereafter  e  is  dropped  whenever  possible.  The  propagation  constant 
is 


•y  =  Q  +  jp 


or,  in  general, 

y(w)  =  a{w)  + 

as  y  is  a  function  of  w  alone  (with  angle  9  and  all  other  parameters  fixed). 
Of  course,  phase  velocity  and  group  velocity  (in  the  "q"  direction)  are. 


respectively. 
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and 

V  =  ^ 

g  ap  ■ 

A  basic  set  of  simultaneous,  linear  differential  equations  can  be 
derived  which  applies  to  a  mixture  of  three  components  and  includes  the 
effects  of  pressure  and  collisions.  The  eigenvalue  solution  of  the  determi¬ 
nant  (of  the  homogeneous  set)  can  be  very  complicated,  so  usually  neutrals 
are  disregarded.  Further  simplifications  omit  pressure  gradients;  further, 
collisions  are  ignored.  As  consideration  of  an  arbitrary  0  gets  difficult, 
the  directions  0^0*  and  90*  are  chosen.  Then,  too,  sometimes  the  dis¬ 
placement  current  term  in  Maxwell's  equations  is  ignored,  only  low  frequencies 
are  considered,  etc. 

At  least  three  considerations  force  these  simplifications.  One  is  the 
tremendous  amount  of  algebra  otherwise  encountered.  A  second  is  the 
obvious  /  correctness  of  the  approximations  in  special  cases.  A  third  is  the 
impossibility  of  exact  calculation  without  a  computer  if  the  approximations 
aren't  made. 

This  paper  shows:  (1)  the  derivation  of  the  basic  set  of  equations 
and  (2)  solutions  of  the  equations  in  certain  cases  of  increasing  complexity. 
Lastly,  there  is  a  brief  discussion  of  research  requiring  knowledge  of 
plasma  propagation  phenomena. 
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II.  DERIVATION  OF  BASIC  EQUATIONS 


We  begin  with  the  macroscopic  approach  and  Boltzmann's  equation 


+  V  •  V,d  +  —  •  V„d  =  Ir 


at 


m 


at 


coll. 


where  "d"  is  the  particle  distribution  function  in  phase  space, 


d  =  f(x,  y,  z,  V  ,  V  ,  V  ,  t ).  Spitzer  shows  how  the  equation  of  momentum 
X  yr  z 

transfer  and  the  equation  of  mass  continuity  are  derived  from  this  starting 
point.  Chapman  and  Cowling  show  that  for  a  gas  of  any  number  of  compo¬ 
nents,  an  equation  of  momentum  transfer  can  be  written  for  each  component, 

.  ,  ,  .th  3 

say  for  the  j  component. 


m. 

J*- 


av. 


=  p.  (E  +  V.  X  B)  -  V  •  a.' 
>  )  1 


-  fm.  i 

^  h 


P. 


jh 


(1) 


where 


V. 


n^m,  =  mass  density 
avg.  velocity  at  R  of  gas 


P: 


n.Z  =  charge  density  (Z  =  1  for  singly- ionized  gas) 
J  ® 


D  -  stress  tensor  =  VP  in  isotropic,  inviscid  case  (P.  =  kinetic 
gas  pressure)  ^  ^ 


I 
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P 


jh 


] 

m. 

] 


gradient  of  gravitational  potential 

momentum  transferred  to  particles  by  "b"  particles, 
particles/unit  volume 
mass/particle 
electronic  charge 


Next,  this  equation  is  specialized  to  a  mixture  of  ions,  electrons  and 
neutral  particles  v^here 


1. 


so  that 


and 


and 


n. 

1 


n^  =  n  =  density 


particles 


m 


at  a  point  (in  a  differential 


volume) 


2. 

n 

n 

3. 

e  = 

4. 

=  -  where  (see  reference  4) 


f..  m.m. 


= 


in 


n,. 

m.+  m. 

h 

}  h 

f, 

m  m. 

le 

e  1 

n 

m  +  m. 

e  1 

f 

m.m 

in 

i  n 

n 

m.+  m 

n 

1  n 

f 

m  m 

en 

e  n 

n 

m  +  m 

n 

e  n 

m 


nn. 


m. 


in 

n  Z 
n 


=  f- 


in  n 


m 


=  f 


en  m 
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if  we  assume  m.  (ion  mass)  is  equal  to  m  (neutral  mass),  and 


m 

1  +  a  unity, 
n 


(Here  the  "f's"  are  effective  collision  frequencies.) 


For  ion,  electrons  and  neutrals,  equation  (1)  becomes 


Ions: 


-  - 

nm.-^+nm  f.  (V. -V  )  +  nm.  f.'  (^.-7  )  =  nq(E+V.  xkB  )-Vp. 
1  9t  e  le  1  e  i  in  i  n  i  o 


(2) 


Electrons; 

av 


nm 


—  *  -  -  »HE  +  V,  .kB^)  -  ,3) 


Neutrals: 


n  m 
n  n 


aVn  ^ 

-~  +  nm.f.'  (V  -V.)+nm  f  (V  -  V  )  =  >  Vp 
at  1  in  n  1  e  en  n  e  '^n 


(4) 


Hoping  we  do  not  invalidate  grandiose  ambitions,  we  have  dropped  the  gravity 

1-^  rm 

and  (V  •  V)  V  terms,  in  order  to  avoid  non-linearities  and  excessive  algebraic 
complications,  as  is  commonly  done. 

In  these  equations,  it  is  necessary  to  note  with  great  care  that: 

!•  E  =  E^  +  5E  =  6E  =  steady-state  field  plus  the  small  perturbation 
(E^  is  assumed  zero,  fiE  varying  as  e  ^^). 
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2.  B  =  -I-  £B  (kB^  »  dB,  so  6B  is  dropped  from  equations  (2), 

(3)  and  (4).) 

3.  n  =  N  .  5n  .  =  undisturbed  charge  density  plus  perturbation 

GO  ^  \Q  G I  X 

(N  =  N.  =  N  hereafter,  but  6n  dn.  necessarily), 
eo  lo  o  e  1 

4.  n  =  N  +  6n  (N  »  6n  ;  N  =  N  hereafter). 

n  n  n  n  n  n 

5.  P  .  =  P  .  +  6p  .  =  P  +  6p 

e,i  eo,  lo  e,i  o  e,i 

=  steady  state  pressure  plus  perturbation 

(P  =  N  kT,  assuming  T  =  T.  =  T  and,  also  VP  or  VP  ,  i.e., 
o  o  e  1  n  on 

space  derivatives  of  unperturbed  pressures,  are  zero  in  the 
homogeneous  layers.) 

6.  V  ,  =  V  +  6V  .  =  6V  .  (The  gases  are  at 

e,i,n  eo,  to,  no  e,  i,  n  e,  i,  n  ° 

rest  and  we  have  dropped  the  "6"  from  equations  (2),  (3)  and  (4). ) 

7.  P=P  +6p=P+dp  with 

n  no  “^n  “^n 

(P  =  N  kT  and  VP  =  0) 
n 

Continuing,  we  use  equations  (2),  (3),  (4)  and  the  mass  continuity 
L  .th 

equation  for  the  j  component. 


-^+  N.  V  •  V.  =  0 
at  j  j 


along  with  the  following  equation  of  state  for  a  perfect  gas,  derivable  assuming 
a  monomolecular  gas  with  "I"  degrees  of  freedom -and  adiabatic  compression. 


8 


dn.  =  8p. 
)  J 


N. 

-1£ 

P. 

)o 


(6) 


the  use  of  which  with  equation  (5)  can  eliminate  the  perturbations  in  number 
density  in  favor  of  perturbations  in  pressure,  for,  substituting  equation  (6) 
in  equation  (5),  we  get 


or 


(7) 


dropping  the  6  notation. 

We  now  proceed  to  derive  a  "generalized  Ohm's  Law"  for  the  ternary 
mixture.  The  pattern  is  fairly  familiar  but  certain  steps  need  explicit  presen¬ 
tation.  Multiply  equation  (2)  by  +  q  and  divide  by  m^;  multiply  equation  (3) 
by  -q  and  divide  by  m^.  Add  the  results  to  get 


Jl_ 

at 


qnV.  -  qnV 
i  e 


m 

+  qnf  {l  +  -^)(V.  -V  )+  qn 
le  m.  1  e 


f;  (V  -  V  )  -  f  {V  -  V  ) 
in  i  n  en  e  n 


l. 

2 

r  1 

J 

1  n 

2 

1 

r  V.  V  _ 

L 

r^Pe 

<1 

-J 

=  q  n 

E  +  q  n 

^  +  * 

X  kB  +  q 

Lm, 

m 

e  -1 

m.  m 

L  i  e-i 

o  ^ 

L""e 

mi  J 

Letting  J  =  qn  (V.  -  V  )  =  N  q  (V,  -  V  ) 
i  e  o  i  e 


9 


m.  m  m 
1  e  e 

tn 

e  1 


and  multiplying  by 


we  have, 


m 


N  q 
o 


^  i  +  f.  J  +  qN 
2  1  8t  le  ^  o 


f.'  (V.  -  V  )-f  (V  -V  )1\ 

in  i  n  en  e  n  J 


r-m  V.  +  m.^  ,  f-  nti 

f-s-j - L_s.  +“rvp  -  — Vp. 

m.m  o  N  q  m.  “^i 

L  ie-l  o  L  iJ 


=  E  +  m. 


Examining  the  second  term  on  the  right,  we  have 


m. 


m  (V.  +  V  -  V  )  +  m.  (V  +  V.  -  V.) 

e'i  e  e  i  e  i  i 


X  kB 


m. 


(V  -  V.)(m.  -  m  )+m.V.  +  m  V 
'e  6  ii  6C 


X  kB 


nm.V,  +  n  m  V 

(V  -v.)+  — 

e  1  n(m.  +  m  ) 

i  e' 


X  kB 


and. 


(8) 


defining 


10 


n.tn.V.  +  n  m  V  nm.V.  +  nm  V 
y  -  ill  e  e  e  ^  ^  *  * 

p  n(ni.  +  m  ) 

i  e 


nm. 


we  have  this  term  equal  to 


(V  -V.)+V 
e  1  p 


X  k  B  or 
o 


■  +  V  lx  kB 

L  N  q  pj  c 


(9) 


Now  we  examine  the  bracketed  terms  in  the  left-hand  brace: 


f!  (V.  -  V  )-f  (V  -  V) 
in  i  n  en  e  n 


in 


m  m 

+  — (V  -V  )  +  —(V.-V.) 

I  n  m.  e  e  m.  i  i 

i  i 


which  becomes 


-f 


en 


m 

V  -  V  +  —  (V  -  V  )  +  V.  -  V. 
e  nm.ee  i  i 
1 


r  nti  m 

=  £.’  (1-— )V. -V  +  — 2-  (V  -V  +V.) 

in  m.  i  n  m.  e  e  i 

^  i  1  -I 


-  f 


en 


m  m 

(1 - -)Y  -V  +— (V)  +  V.  -V. 


m. 


(f  -  f ) 

en  in 


L 

i 

i 

1- 

e 

V  ■  V 

X  ^7^ 

£  + 
en 

m 

m. 

L  i 

V  V  . 

e  I 

T  V 

n 

+ 

—  V 
m.  in 

i 

1  i 


(V.  -  V  ) 
1  e 


£  -£.' 
en  in 


n  m  V  +  m.  V. 
e  e  i  i 


m. 

i 


-  V 


m 

£  +  ~  f ; 

en  m.  in 
i 


(V.-V^) 


11 


where,  letting 


(f  -£.')  =  f,  . 

en  in  1 


we  have 


-f. 


r  1  r  n 

V  -  V  +  £  +  ~  f! 

p  n  en  m.  in 

L  J  L  1  J 


(V.-VJ 


(10) 


We  now  substitute  expressions  (9)  and  (10)  into  equation  (8),  with 


m 

£.  +  f  +  — 2.  £.  )  =  £, 

le  en  m.  in  y  3 
i  ' 


to  get 


m 
_ e 

q^N 


If  -3^ 


m  £.  ■? 

(V  -V  )  +  ;:^  X  kB  =  E  + V  X  kB 
q  P  n  N^q  o  p  o 


and,  since 


«  = 
e  m 


the  cyclotron  resonance  £requency  £or  electrons,  we  obtain  £inally 


1 


2 

«  u 
o  p 


3J  ——  — 

— +£_J+wJxk 
9t  3  e 


-  -  ®  ■ 

E  +  V  xkB  + (V  - 

p  o  u)  p  n 


N  q 
o 


m 

Vp - 2.  Vp. 

“^e  m.  *^1 

i 


(U) 


12 


5, 6 

Equation  (11)  is  the  general  expression  sought.  We  have  introduced 
the  electron  plasma  frequency 


u 


2  No'l 


p  m  ( 
e  o 


m 


The  equation  is  in  the  RMKS  units.  No  approximations  other  than  l±  —  =  1  and 

i 

m.  =  m  have  been  made. 

1  n 

We  now  proceed  to  find  the  equation  of  motion  for  (1)  the  charged 
component  and  (2)  the  neutral  component.  Returning  to  equation  (4)  on  page  6, 
we  have 

av 


0  + 
*^m  8t 
n 


+  Prr.  ‘ 

m ,  in  n  1  m  n  e 
1  e 


=  -  Vp. 


(4') 


and,  examining  the  bracketed  term  as  before,  making  similar  approxima¬ 
tions,  we  have 


m 


f.'  r  V.  -  V  + — 

.mil  n  m. 

i  i 


m 


p  t;  I  V.-V  + - (V. -V.)+— ^(V  -V  ) 

~  ~  1  1  m.  e  e 

1 


m 

+  p  f  (^  -  7  +  V.  -  ~  (V  -  V  ) 

menn  e  i  im.  e  e 

e  i 


p  [f  -  f.'  ]  (V.  -  V  )  -  p 
m  '■  en  in'  i  e  m, 
e  i 


m 

f!  +  ~  f 
in  m.  en 

i 


m 

V.  +  ~V  -V 
i  m.  e  n 

i 


which  may  be  written,  recalling 


f  -  f.'  =  f, 

en  in  t 


L3 


N  q  (V.  -  V  )  =  J 
o  x  e 


^m.  V.  +  V 
11  e  e 


N  m. 
o  1 


and  defining 


m.^.  +  m  y 

^  =_LJ - !_£.  a 

P  m. 


^2  = 


f  + 
in 


£  \ 
m.  en  J 
1  ^ 


as 


tn 


—  f,  J  -  N  m  £  (V  -  V  ) 

q  I  o  i  2  p  n 


which  is  also 


B  £  -*  _* 

-2-^  J  +  N  m.f^  (V  -V  ). 
u  o  1  2  n  p 


(12) 


Substituting  equation  (12)  into  equation  (4')  we  obtain 


Nm 


av 

_ r 

i  at 


B  £ 

-2-i-  J  -  N  m,£,  (V  -  V  )  -  Vp 
w  o  1  2  n  p  n 

e  ^ 


(13) 


To  obtain  an  equation  o£  motion  for  the  charged  component  o£  the 
gas,  we  add  equations  (2)  and  (3)  to  get 


av 


9V. 


Prr.  TT  +  ^  +  Pm  ■  'V„)+  P,^  £^^( )  =  J  X  B 

m  at  m.  at  m,  in  i  n  m  en  e  n  o 


-  Vp,  -  Vp 
1  e 


14 


or 


8^ 

8t 

P 


(V  -  V.)  +  p  f  (V  -V  ) 
in  n  i  m  en  n  e 
I-  i  c 


J  X  kB  -  Vp.  -  Vp 
o  1  e 


The  bracketed  terms  have  been  shown  equal  to  expression  (12);  so,  on 
substitution,  we  have 


8V  B  ^1  ^ 

N  m.  =  -2-1- J  -  N  m.f  (V  -  V  )  +  J  x  kB  -  Vp  (13') 

oidt  u  oi2  pn  o  ‘^p 


Next,  from  the  equation  of  mass  continuity,  in  the  form  of  equation  (7), 
we  can  write  for  the  electron  gas 


p  \i*zj  at  e  • 


This  can  be  written  as 
1  /  I  \  9Pe 


1+2  J  8t 


+  V 


V  -  rr“ 

p  N^q 


=  0 


(14) 


since 


V  - 
p  N^q 


m  V  +  m.V,  N  q  (V.  -  V  ) 
ee  li  o  i  e 


m. 

1 


N  q 
o 


m 

—2  V  +  V.  -  V.  +  V 
m.  e  I  I  e 
1 


m  . 

1  +  — 2  )  V  a  V  . 

e  e 
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For  the  charged  gas,  we  have  by  adding  equation  (7)  for  the  electron 
and  ion  gases  together 


8p 

_LE  +  V 

at 


(15) 


since 


and 


P  =  N  k  T 
o  o  e,  1 


2P  =  (N  +  N  )  k  T  . 
o  eo  lo  e ,  1 


-  7 

IV  -  =  2~V  +  2V.  -V.  +  V 

pNq  m.  e  lie 

o  i 


.  m  \ 


V  +  V. 
e  1 


=  V  +  V. 
e  1 


For  the  neutral  gas,  we  have 


,  /  f  X  9p  _ 

^  ^  —S.  +  V  •  V  =0 

PVf  +2^  at  n 


(16) 


In  conclusion,  we  write  Maxwell's  equations  as 


V  X  (V  X  E)  =  -  jwft  J  +  E 


(17) 


16 


and,  collecting  equations  (11),  (13),  (13*),  (14),  (15),  (16)  and  (17)  together, 
letting  =  m,  we  have 
6hm'''.a  Law: 


c  w 
o  p 


(J,.  t  tj) j  t  X iT]  =  E  +  5^  X *  "iv 


Equation  of  Motion  for  the  Charged  Component: 

jwN  m  V  =  J  -  N  m£,  (  V  -V  U  J  xkB_  -  Vp_  (13') 

■^o  poj  o2\^p  ny 


o  p 


Equation  of  Motion  for  the  Neutral  Conniponent: 

^  ®  ^  ^  ^  X 

juNmV  =  -  J  -Ntnf  (  V  -V  )  -  Vp 

J  n  u  o2\^npy  n 


(13) 


Continuity  Equation  for  the  Electron  Gas: 


- £.  +  V  •  f  V  -  — —  )  =  0 

P  k  ^  I  ''p  N  q  ' 

o  N  O 


(14') 


Continuity  Equation  for  the  Charged  Gas: 


2P  k 
o 


y  .  /  V - - - ^  =  0 

V  P  2N  q; 


(15) 


Continuity  Equation  for  the  Neutral  Gas: 


]  wp 


n 


+  V  V  =  0 
P  k  n 


(11) 


(16) 
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Maxwell's  Equation*  (Wave  Equation): 

Vx^VxE^  =  -  jufi  J  +  fi  E  (17) 


The  statement  of  the  last  seven  equations  concludes  this  section. 

Derivatives  with  respect  to  time  have  been  taken.  All  p^'s  (pressures) 

—  .  1  /  f 

are  perturbations  only;  kB^  is  the  fixed  field define  —  =  < - 


\l+  2 


) 


and 


1 


J  +2 
n 


m 


We  have  supposed  that  Vp  »  -  ^p.-  These  seven 

e  m.  1 
1 


equations  have  seven  unknowns,  which  are  written  J,  V  ,  V  ,  p  ,  p., 

p  n  e  I 

p  I  E.  These  seven  equations  are  used  in  the  next  section.  The  reader 
n 

would  do  well  at  this  point  to  refer  to  Watanabe's  thorough  work  on  this 

,  .  5,6 

analysis. 


III.  SPECIAL  CASE  SOLUTIONS 

A.  The  MHD  Approximations  and  a  Simple  Case; 

As  an  illustration  of  the  use  of  the  basic  seven  equations  and 

to  show  in  a  most  simple  way  the  method  to  be  used  later,  we  make  what 

1  2 

are  called  "magneto -hydrodynamic  (MHD)  approximations".  ’  We 
consider  only  a  fully-ionized  gas,  with  no  neutrals,  keep  only  one  pressure 
term,  set  all  collision  frequencies  equal  to  zero,  disregard  displacement 
current  to  get  the  following  set: 


18 


E  =  -  V  X  kB 

P  o 

(18) 

iu.'4s  iBlr,  V  =JxkB  -Vp 
•op  op 

(19) 

VxVxE=-jwMJ 

(20) 

■— —  p  V  •  '7  =0 

2  P  k  ^p  p 

(21) 

We  first  obtain  from  this  set  two  equations  in  E  and  by  solving  equation 
(21)  for  and  substituting  into  equations  (19) >  and  solving  equation  (20)  for 
J  to  substitute  into  eq.  (19)>  We  get 


E  =  -  V  X  kB 
P  o 

^  B 

juN  mV  =--r^  VxVxE 

o  P  1 


^  2  P  k 

X  k  +  ■  V(V  •  V  ) 

}  P 


We  now  let 


E  = 


i  E  -)■  i  E  +  kE 

X  y  1 


y  X  sin  $  -  y  z  cos  9 


and 


V  = 
P 


i  V  +  j  V  +  kV 
X  y  z 


-  y  X  sin  9  -  y  z  cos  9 


(22) 

(23) 


substituting  into  equations  (22)  and  (23),  separate  each  equation  into  i,  j,  k 
components,  and  we  get  six  equations: 


E  +  B  V  =  0 
X  o  y 


E  -  B  V  =  0 
y  ox 
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E  =  0 
z 


B  ,  2P  k  .  ,  _ 
...  o  o  f  Z  £ 

iwN  mV  - V  E  -  ' "  (  y  o, 
■'  o  X  juft  y  ]u  I 


,  2P  k  , 

V - rr —  a.o.  y  ^  =  0 

X  ]cp 


®o  r  2  2  2 

iuN  mV  +T~‘  y  o..  E  -y  o, a,E 
o  y  ]WM  I  2  X  I  2  z 


■^r 

2P  k 

jwN  m.  V  — y  o.a.  V  +  y  a.  V 
■’oiz  juL'12x'2 


=  0 


where  we  have  put  sin  6  =  a  and  cos  6  =  a  ,  these  not  to  be  confused  with 

^  M 

attenuation  factors. 

We  know  that  for  a  solution  to  exist  the  determinant  of  the  co¬ 
efficients  must  be  identically  zero.  This  is 


B 


1 

0 

0 

0 

o  2  2 


■a_  y 
jwfi  2 


0 

I 

0 

B  y2 
o 

0 


0 

0 

1 

0 


0 

■B 


2P  k 

jwN  - - o  y 

o  jw  1  ' 


®o  2 


+B 

c 

0 

0 

0 


juN  m 
o 


0 
0 
0 

2 

ao  y 
ju  12 

0 


2P  k 
o 


2P  k 
o 


any 
jw  12 


2P  k  ,  , 
o  2  2 
— : a,  y 

jw  2 


j  wN  m 
o 


Most  conveniently,  an  expansion  leads  to  the  following  factors: 

B^  ,2P_k  ,  ,  ,  ,2P  k 


.  _ o  ,  .2 

jwN  m  -  *: - ( a_y) 

'  o  jwfi  2' 


JWfJ 


/  o  2  2  \  /  o  2  \ 

\“Tr’“2'»'  -jwN^mj 


(24) 
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Thus  we  see  three  modes  of  propagation,  for  an  arbitrary  angle,  one 
mode  for  0  =  90*  and  two  for  0  =*0*.  For  0  =  0*,  o  =  0  and  a  =  I,  and 

i  M 

we  have 


P  N  m  _ 

1 

p2P  k 

jw  j 

o 

\f  — 

o 

2P  k 

L-  O  -1 

1  ^ 

P\- 

N  m 
I-  o 

1 

2 


N  m  ^ 


JO) 


O 


o 

N  m 
■-  o 


1 

2 


where  y  represents  a  sound  type  or  electro-acoustic  wave,  y  the  well- 

a 

Z 

known  Alfven  mode,  where  /p  ,  the  magnetic  pressure,  corresponds  to 
2  kinetic  gas  pressure. 

We  see  that  for  the  arbitrary  angle,  y^  and  V  are  found  by  simply 


substituting  for  B  ,  B  cos  .0« 
o  o 

For  0  =  90°,  a  =  0  and  o ,  =  I  and  we  have 

C»  a 


jw 


N  m 
o 


2P  k  +  B^/iu 
o  o 


ph3  = 


2P  k  +  B„/n 


N  m 
o 


and  appropriate  changes  are  made  for  an  arbitrary  angle.  The  kinetic  and 
magnetic  pressures  both  clearly  affect  the  third  mode. 

2 

In  many  physical  situations,  one  of  the  terms  B  or  2  P  k  is  much 

o  o 

greater  than  the  other  and  expressions  simplify. 


21 


B.  Additional  Cases  Factorable: 

As  long  as  the  6x6  determinant  obtained  from  the  E  and 
equation  is  factorable  (and  this  sometimes  requires  some  experimentation  to 
find  the  right  expansion),  the  y's  for  the  various  modes  can  be  solved  for 
without  too  much  trouble.  Suppose  we  have  a  set  of  equations  thus: 


Jl±. 


2 

c  w 
o  p 


E  +  V  X  B  k 
P  o 


iu  N  m 
o 


V  = 
P 


J  X  k  B  -  Vp 


2P  k 
0 


p  +  V  •  V  =  0 
P  P 


8J 

Here  we  have  included  the  displacement  current  and  8t  ,  but  have 

ignored  all  other  terms  and  collisions  in  the  basic  set  of  seven. 

The  two  equations  for  E  and  V  are: 

P 


i<j 


€ 


I 

iwu 


("vxVxE-  E^]  =  E+  V  xk 


ju  N  m  V 
o  p 


B 

_ o_ 


V  x  V  X  E 


2 

u  u  c 

o 


E 


-  ^  2P  k 

X  k  +  V(V  •  V  ) 

P  jw 


and 
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In  this  case,  after  the  6x6  determinant  is  found  and  expanded 


we  get: 


1.  For  any  angle  9:  (c  =  velocity  of  light  in  vacuum) 

jwN^m^ 


,  ^  ,2P  k  - 

juf  B  -  7—^  —  (o  y) 

J  o  o  jufi  ]  w  '2 


2Pok  2  N 

+  jwN^m(— T~y  -jV)N^mj 


o  V  ju 


2  2,, 
u  c  2 

y 


\  OJ  w 


p  p 


)\ 

' 

(,  V 

L\  / 

P 

,  ,-2  2, 

Bo  ,  i2  .  f .  w  c  2 

T**"  (a.  y)  -  iwN  m  (  1  -  —r - -  y 

Ju^l  '  2  "  o  V  2  2  ' 

P  P 


2  2 


_  2  r  u  c  , 

-  3  w  €  B  1 - -  -  — r  (o 

■'  o  o  2  2  1 

L  CO  U 


Y) 


(25) 


P  P 


which  for  fl  =  0*  becomes 


/  2  B^  2\ 

(jw€  B  -r —  Y  |+j«N 


2  .2  ,_2,2Pk 


,,  M  c _  2 

U  CO 

p  p 


o  2  .  „ 

- y  -  ]wN  m 

jw  o 


=  0  (26) 


which  suggests  how  the  Alfven  mode  is  affected  and  shows  again  the  common 
acoustic  mode. 

2.  For  6  =  90",  we  have 

y  )+(jcoI 

o  o  jw^  y  V  °  J  “ 


->  n“  >  ✓  2P  k  2  2 

jw  -  ry^y  j+f  jcoN^m - r;-; —  y  )( 1 


±L  ^  2\- 

2  ■  2  ) 
w  w  ''•J 


P  P 


2  2  , 
,  w  c  2 

1---  —  y 

CO  w 

P  P 


0. 


(27) 
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Explicit  expressions  for  the  y's  have  not  been  solved  for  as 
this  and  the  next  case  are  regarded  only  as  "pattern-setters"  for  preparation 
on  the  more  complex  cases. 

As  the  next  case,  suppose  that  the  given  set  of  equation  is 

the  following: 


I 


2 

t  <ti  ^ 

o  p 


iuj  J  X  k 
•’  e 


=  E  +  V  X  k  B 
P  o 


iwN  mV  =  JxBk  -  Vp 
op  o  p 


V  X  V  X  E  =  -  jupL  J  +  W^/X€^E 


jwP 


rrs-t-  +  V  •  V  ,=  0. 
2P  k  p 

o 


This  set  has  been  solved  in  the  same  manner  as  before  and  it  is  found  that 
the  determinantal  expansion  does  not  factor.  It  has  thus  been  tentatively 
concluded,  and  it  is  most  probable,  that  no  factors  exist  in  the  complex 
cases  of  physical  interest.  The  next  case  treated  moves  us  farther  into 
this  realm. 

C.  The  FuUy-Ionized  Case 

We  now  consider  a  fully-ionized  plasma,  and  use  the  basic 
set  of  equations  with  collision  terms  set  to  zero  (tenuous  case).  Our 


basic  set  then  is; 
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I 


2 

(  u 
o  p 


ju j  -t-  X  k 


=  E+  V  xkB  +  „ 

p  o  N^q 


^Pe 

(28) 


juN  mV  =  J  X  kB  -Vp 
op  o  p 


(29) 


V  X  V  X  E  =  -  jup  J  +  E 


(30) 


^  p  +  V  .  (^v  =  0 

P  k  Pe  p  N  q; 


(31) 


lU 


p  +  V  •  (  V 

2P  k  ^p  I  p  2N 

o  \  o 


(32) 


The  two  equations  for  E  and  are: 


_  w  .. 

I  -iL. 

g  + 

r  i  1 

2 

u  J 

Lfn  u  J 

V  X  V  X  2  +T 


]wp 


u 


2 

U  *  w  J 

o  p 


X  V  X  E^  X  k 


+  jw 


r- 

e 
2 

u  -I 
P 


k 


Exk+V  xkB 

p  o  j  u 


(’  •  %) 


«  P  k  ^ 
o  o 

m2  2 

N  q 
o 


V  (V  •  E)  =  0. 


(33) 


and 
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B 

o 


(VxVxE)xk  + 


iuc  B 
o  o 


Exk+juN  mV 
o  p 


r  n  p  *  * 

2  P  k  V  (V  •  V  )  - 

o  j  '  p  L  I 


c  P  k 
o 


V  (V  •  E)  =  0.  (34) 


In  the  same  way  as  previously,  we  put  in  the  exponential  forms 

for  £  and  and  find  the  six  homogeneous,  simultaneous,  linear  equations. 

Next  the  6x6  determinant  of  the  co-efficients  is  set  equal  to  zero.  This 

determinant  is  shown  on  page  26  ( 

The  solution  of  the  eigenvalue  problem  is  an  interesting 

exercise.  First,  for  6  =  0*,  the  determinant  is  factorable  into  two  quadratic 
2 

equations  in  y  and  the  four  modes  are  as  follows: 


N  m 
o 

2P  k 
o 


2 

+  u> 


N  m 
o 


P  k 
o 


and 


2 

■^3 


and 
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I  + 


(j 


2  2 
c  u  N  m  <d 
o  p 


74  = 


B 


w^uN  m 
o 


i  “  V 


It  will  be  noticed  that  Y  ,  delated  to  the  original  acoustic  wave  but 

l»  2 

modified,  and  Yo  a  related  to  the  original  Alfvdn  modes. 

J  I  ^ 

As  a  check  on  these  results,  let  us  suppose  that  =  0. 
and  becomes 

'3  U 


(j 


u 


u 


1  - 


b> 


U 


'3,4 


w 


Therefore 


]W 


2 

w  I 

«  -> 


'3,4 


and 


ph 


3,4 


I  - 


w 


w 


This  is  a  familiar  result. 


1,2 


is 


Then 
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V 

Suppote  further  that  the  ratios  and 
where  is  the  simple  Alfvdn  velocity,  and  we  obtain 


V 

— ^  are  small, 
u 


80  that 


and 


w 


1 


1 
n  2 


2  w 


These  two  results  are  also  familiar. 


1,2 


Ph, 


is  the  phase 


velocity  of  the  "extraordinary"  wave  and  that  of  the  "ordinary"  wave. 
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Figures  2  and  3  ehow  the  variation  in  and  y^  with  w,  for  the  complete 
expressions  on  page  28.  Of  course. 


<1B, 


u. 


i  m. 
1 


and 


N  q 

w  =  —2 — 
p.  ni.c 
*^1  1  o 


When  y  is  positive,  the  wave  does  not  propagate,  but  is  evanescent.  The 
y^  mode  iS  essentially  connected  with  the  ions  and  the  y^  mode  with  the 
electrons. 

We  shall  not  consider  2  acoustic  modes.  We  refer 

7 

the  interested  reader  to  Pai's  paper. 

Next,  we  may  examine  the  case  0  =  90*.  When  we  let 

V> 

a  =  0  and  a  =  1  in  A  and  solve,  we  find  one  mode  immediately: 


2 

w 

2  _  Je 

^1  "  2 


1  - 


h) 


[‘-4] 

L  (*)  J 


U) 


b) 


which  is  the  electromagnetic  mode  found  before  the  0=0*  case,  page  Z7. 
However,  this  simple  factor  is  the  only  one, possible  and  we  are  left  with 


a  cubic  equation  in  y  . 
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This  case  as  well  as  the  0  =  0*  case  is  discussed  extensively 
8 

by  Tanenbaum.  However,  the  case  for  arbitrary  0  is  a  very  difficult 
problem  and  has  not  been  solved,  and  it  appears  to  be  a  computer  problem 
for  special  cases. 

This  problem  has  up  to  now  disregarded  electron  and  ion 
collisions.  Usually  this  is  taken  as  considering  infinite  conductivity. 
However,  if  f^^  is  included  in  the  basic  set,  equation  (28)  on  page  24  is 
written  as  ' 


1 


2 

c  u 
o  p 


(ju  +f)J  +  w  Jxk 
3  e 


^  ^  ""Pe 

=  E  +  V  X  kB  + 

p  o  N  q 


m 

where,  instead  of  f,  being  (f  +  f  +  f.  )  as  in  the  basic  set,  it  is  only 
3  le  en  m,  in 

1 

f.  .  This  insertion  changes  the  6x6  determinant  very  little,  but  enough 

to  introduce  complex  y  's  so  that  attenuation  of  the  modes  is  caused.  It  is 

known  that  at  very  -  low  frequencies,  MHD  effects  hold  true  and  at  very 

16 

high  frequencies,  the  Appleton-Hartree  equations  can  now  be  derived. 
Computer  calculation  of  attenuation  and  phase  factors  for  cases  of  special 
interest  is  possible,  as  before. 


D.  The  Partially-Ionized  Case: 

In  the  ionosphere  and  in  most  man-made  plasmas,  it  cannot 
be  assumed  that  full-ionizatioti  exists.  To  study  this  case,  one  must  employ 
each  of  the  basic  seven  equations  on  pages  16  and  17.  We  use  the  same 


32 


method  ae  before,  i.e. ,  solve  for  two  equations  in  E  and  V  ,  and  find  the 

P 

eigenvalues  of  the  6x6  determinant  obtained  after  substitution  of  the 
exponential  forms.  . 

This  is  a  long  and  demanding  project  with  most  excellent 

chances  for  algel^raic  errors.  We  will  present  the  two  equations  in  £  and 

V  ,  which  become  equations  (33)  and  (34)  if  N,  the  neutral  number  density, 

P 

is  set  to  zero  and  all  collision  terms  -  f^,  f^,  f^  -  are  also  zero,  (as  would 

be  expected  if  the  equations  are  correct).  These  equations  for  the  partially- 

ionized  case,  including  the  effects  of  collision  and  pressure  terms,  are  on 

pages  33-34,  The  6x6  determinant  is  on  page  35  (A^).  Three  independent 

checks  have  been  made  of  these  very  complex  expressions. 

Once  again  we  are  presented  with  the  problem  of  solution  for 

0  =  0”,  0  =  90*,  and  arbitrary  0.  First,  it  can  be  seen,  after  much  algebra, 

that  for  0  =  0*,  the  determinant  factors  into  two  polynomials,  one  a  quad- 
2  2 

ratic  in  y  and,  another,  a  cubic  in  y  .  There  are  thus  five  modes  of 
propagation. 

2 

Secondly,  for  0  =  90*,  the  factors  are  a  first  degree  y 
2 

polynomial  and  a  quartic  in  y  (  and  a  very  complicated  one).  Again  we  have 
five  modes. 

Tanenbaum  has  done  an  excellent  analysis  of  these  two  special 

10 

cases  and  has  shown  qualitatively  some  very  interesting  results.  For 
0  =  0*;  the  variation  of  the  roots  of  the  quadratic  are  as  shown  in  Fig.  4, 
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(A) 


]  V  [  v.(  vvp)  )  ]  ^-  [  JhTmV.u  ]  ^  {v.[(V.(».E).  k  ]} 


[i-k/-x!i,^(P»q/)]  [^.QgC  )]  [  L-y^fi  M*Ry^)  ]  0  [  -^^agC  R/)] 
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The  following  statements  apply  to  Fig.  4; 

1.  is  the  phase  velocity  of  the  usual  left - 
circularly-polarized  wave  associated  with  the  ions,  that  of  the  right- 
circularly-polarized  wave  associated  with  the  electrons.  Resonances  occur 
at  the  cyclotron  frequencies. 

2.  For  very  low  frequencies  (MHD),  the  only 
propagating  mode  is  the  Alfvdn  mode  associated  with  the  entire  gas,  and 


V 

a 


_ o_ 

_ H _ 

Total  Mass  Density 


1. 

2 


3.  As  the  frequency  reaches  the  neutrals  cease 
to  take  part  and  we  have  an  Alfvdn  mode  in  the  charged  gas  only 


_ H _ 

Charged  Mass  Density 


1 


4.  As  the  frequency  continues  to  increase,  we  have  the 
Alfvdn  modes  ceasing  to  be  effective  and  we  pass  to  the  magneto-ionic  limits 
as  in  Figs.  2  and  3.  The  dispersion  relation  becomes  approximately  of 
Appleton-Hartree  form. 

5.  These  curves  are  drawn  without  considering 
appreciable  collisions,  but  when  collisional  damping  effects  are  effective, 
the  waves  can  propagate  with  attenuation  at  all  frequencies.  For  high  values 
of  f^,  f^t  f^  the  resonant  effects  are  greatly  diminished. 
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2 

Next,  for  the  cubic  in  y  ,  the  roots  take  the  form  shown  in 

Fig.  5. 

The  following  statements  apply  to  the  curves: 

1.  The  symbols  U  stand  for  the  velocity  of  sound  in 

the  indicated  gas. 

2.  At  very-low  frequencies,  only  the  acoustic  wave  in 
the  total  gas  in  propagated, 

k  P  +  2kP  A 

n _ o _ 

Total  Mass  Density 


3.  When  f^  is  reached,  coupling  to  the  neutrals 
decreases,  and  we  have  an  acoustic  wave  in  the  charged  gas 

2k  P 


U  = 
P 


1 
T  2 


Charged  Density 


and  another  in  the  neutral  gas. 


U 


n 


N  m 


4.  As  w  is  reached,  we  see  that  coupling  between 
P: 


the  ions  and  electrons  decreases,  so  we  have  simply  U.. 


5.  The  electron  oscillations  are  greatly  attenuated 


unless  u  >  w  ;  the  phase  velocity  can  be  very  high  for  u  3  u  ,  but  decreases 

Pe  Pe 

to  U  as  u  -►  00 . 
e 
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6.  Tanenbaum  states  some  of  the  high-frequency 

results  here  cannot  be  trusted  as  they  might  pertain  to  wave-lengths  (for 

10 

high  enough  w's)  which  are  shorter  than  the  Debye  distance. 

For  the  case  of  6  =  90*,  one  root,  as  has  been  said,  can  be 
quickly  obtained,  just  as  in  the  fully-ionized  case  on  page  29.  The  factor  is 


2  f  LA  -  EH1 
^[LC-JE] 

with  L,  A,  G,  J,  E  drawn  from  the  definitions  in  equations  (B)  and  (A),  page  29. 
This  mode  is  obviously  related  to  the  fully-ionized  mode.  Tanenbaum  affirms 
this  by  showing  Figure  6. 

All  other  waves  propagating  at  0  =  90*  are  coupled  longitudinal 

and  transverse  waves  (the  E  vector  is  at  an  acute  angle  with  propagation 

2 

direction).  This  fourth  order  polynomial  in  y  has  not  been  solved  in  general. 

However,  by  the  use  of  known  limits  and  assuming  continuity  of  the  roots, 

10 

Tanenbaum  has  presented  Figure  7,  although  no  analytic  solution  exists. 

The  following  statements  apply: 

1.  At  low  frequencies  only  a  coupled  acoustic  and  MHD 

wave  can  propagate.  This  was  shown  in  the  simplest  case  treated  on  page  20, 

where  , 

V  . 
phase 


2Pk  +  B  //i 
o  o 


TOTAL  DENSITY 

and  an  equivalent  expression  can  be  written  for  the  modified-ion  wave  curve 
in  this  partially-ionized  case  as 
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■o  that  if  the  magnetic  pressure  is  much  greater  than  the  kinetic  gas  pressure, 
the  phase  velocity  is  simply  V  .  the  Alfven  speed  in  the  total  gas, 

2.  As  the  frequency  goes  up,  the  neutrals  are  "left 
behind"  and  we  have  the  appearance  of  the  modified  neutral  wave,  that  has  a 
relatively  large  damping  factor  due  to  neutral-ion  and  neutral-electron 
collisions,  and  the  Alfven  mode  in  the  charged  component  of  gas  V 

3.  As  the  frequency  increases  past  the  effects 
of  the  magnetic  pressure  are  lost  and  we  have  the  acoustic  wave  in  the 
charged  component,  which  next  becomes  an  acoustic  wave  in  the  ion  gas  alone. 

4.  The  remaining  two  modes  are  electron  oscillations  - 
one  longitudinal,  the  other  transverse, 

5.  The  cut-off  frequencies,  or  transition  frequencies, 
may  vary  according  to  the  situation,  but  Tanenbaum  has  written  expressions 
(approximate)  for  them, 

6.  These  curves  apply  when  collision  effects  causing 
damping  are  negligible.  Tanenbaum  discusses  at  length  the  effects  of  damping. 

For  an  arbitrary  angle,  we  have,  of  course,  five  modes  but 

this  situation  has  not  been  solved.  It  is  possible  to  reduce  the  6x6 

2 

to  a  4  X  4  determinant  and  obtain  a  quintic  in  y  .  Numerical  computation  then 
is  demanded,  both  for  the  calculation  of  the  coefficients  of  the  polynomial  and 
the  finding  of  the  roots.  Work  in  this  project  is  now  being  carried  out,  in 
order  to  develop  a  rapid  computer  program  as  well  as  enable  identification 
of  mode  types. 


I 

t 
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IV.  CONCLUSION 

It  is  now  clear  that  algebraic  complications,  even  under  the  assumpr 
tions  of: 

L.  Negligible  gravity  effects 

2.  Equal  electron,  ion,  and  neutral  temperatures 

3.  Negligible  (V  •  V )  V  terms 

4.  Small  amplitude  perturbations 

5.  Negligible  viscosity  effects 

6.  An  approximation  for  effective  collision  frequency 

7.  Spatial  derivatives  of  undisturbed  gas  equal  to  zero 

8.  Adiabatic  compressions 

9.  No  ionization  or  recombination  effects, 

make  this  previously  explained  approach  to  the  problem  of  finding  a  general 
analytic  expression  for  the  various  modal  dispersion  relations  an  impossible 
task,  and  that  only  with  numerical  computation  for  the  case  of  interest  at 
the  particular  moment  can  roots  be  found  and  also  that  one  is  fortunate  if 
the  physical  situation  is  such  that  certain  terms  can  immediately  be  discarded. 
Watanabe  has  done  an  excellent  paper  illustrative  of  this  approximate  or 

intuitive  method  and  includes  tables  for  the  important  physical  parameters  of 

.  5,6 

the  ionosphere. 

Knowledge  of  the  exact  dispersion  relations  will  point  the  way  to 
clarification  of  problems  of  propagation  in  the  ionosphere  under  various  con¬ 
ditions.  Our  problem  has  dealt  with  the  presence  of  a  magnetic  field.  If  the 
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analysis  had  been  performed  without  a  field,  it  would  have  been  possible  to 
uncover  the  fact  that  the  imposition  of  a  fixed  field  in  the  vicinity  of  an  antenna 
on  a  re-entry  space  vehicle  should  help  to  overcome  the  loss  of  radio  contact 
occasioned  by  the  hot,  ionized  plasma  sheath. 

For  the  last  ten  years,  from  various  stations  over  the  world,  great 
amounts  of  recorded  data  have  been  obtained  pertaining  to  micropulsations  of 

the  earth's  magnetic  field.  Amplitudes  range  from  fractions  of  a  gamma 

-5  12 
(=  10  gauss)  to  tens  of  gammas,  with  frequencies  from  0.01  to  10-20  cps. 

The  main  reason  for  studying  the  phenomenon  is  that  the  generation  of  such 

pulsations  is  thought  to  be  related  to  events  in  outer  space.  Dungey  in  1954 

first  suggested  that  geomagnetic  micropulsations  could  be  explained  by  MHD 
13 

oscillations.  A  great  amount  of  work  has  been  done  to  develop  basic 

equations  in  rectangular,  cylindrical,  and  spherical  co-ordinates  to  describe 

14  15 

propagation  in  the  slightly  ionized  gas  of  the  ionosphere.  '  The  effort  has 

been  to  explain  the  observed  amplitudes  and  change  in  character  of  the  fluc- 

l6 

tuations  with  time  of  day,  location,  sun-spot  activity,  etc.  This  research 
is  continuing. 

Another  problem  of  research  interest  requiring  familiarity  with  the 

type  of  analysis  of  this  paper  is  that  of  ionospheric  heating.  The  attenuation 

factor  in  the  propagation  factor  is  of  course  pertinent  here.  Akasofu  claims 

that  the  energy  dissipation  process  of  hydromagnetic  waves  in  the  ionosphere 

and  the  consequent  heating  in  the  F  region  of  the  ionosphere  can  only  be 

17 

negligible,  of  the  order  of  S’C  under  the  most  favorable  conditions. 
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The  study  of  the  propagation  of  signale  from  sudden  disturbances  in 

the  troposphere,  such  as  a  lightning  flash  giving  rise  to  "whistlers",  or  in 

the  ionosphere,  such  as  high-altitude  bomb  burst  or  meteoritic  activity ,, can 

iS 

be  clarified  with  the  use  of  concepts  introduced  in  this  paper. 


>  '< 
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